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C""~~ , Abstract. We discuss the derivation of master equations in the presence of initial correlations 

^~~^ • with the reservoir. In van Hove's limit, the total system behaves as if it started from a factorized 

I initial condition. A proper choice of Nakajima-Zwanzig's projection operator is crucial and the 

►..^ . reservoir should be endowed with the mixing property. 

^H ■ 

tH- ' 1. Introduction 

o : 

\^ . The derivation of master equations has been a central issue for decades in the 

^P I field of quantum physics, both from a practical and fundamental point of view. 

,J^ ' The interests include dissipative dynamics in quantum optics [1,2], understand- 

^~^- ing decoherence [3], the emergence of irreversibility [2,4], and the influence of 

■y i dissipation/decoherence on quantum information. 

C^ ' It is not a trivial problem how to describe the dissipative and irreversible 

^' dynamics of quantum systems, since the Schrodinger equation applies to closed 

• . , systems and describes unitary evolutions. In the standard approach to this issue, 

. 5^ I one puts the system in a large reservoir to allow energy exchange and dissipation. 

/\^ • Then, one applies the Schrodinger equation to the total (closed) system, system 

S+reservoir B, and extracts the dynamics of system S by tracing out the reservoir 
degrees of freedom, 

ps{t) =tVB p{t), (1) 

which exhibits a non-unitary evolution and is described by a master equation. 
Here, p{t) and ps{t) are the density operators of the total system and system S, 
respectively. 

One of the targets of this article is to discuss a hypothesis that is taken for 
granted in most of such approaches: factorized (uncorrelated) initial conditions 

P0 = PS® PB (2) 

are usually assumed in the derivation of master equations without any reasoning 
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about the validity and the necessity of this hypothesis [1-4]. We will discuss in the 
following a derivation of master equations for generic (correlated) initial states. 
The reconsideration of the factorized initial condition (jj} raises another question 
in the application of Nakajima-Zwanzig's projection operator V [2,4-6], defined 
by the mapping 

p{t) -^ Vp{t) = trB{p(t)} ® Ob. (3) 

Which state of reservoir B should be taken as the reference state JIb of the projector 
VI We address these points and show the derivation of master equations in the 
presence of initial correlations. We clarify that the choice of V is crucial and the 
reference state VL-q should be endowed with the mixing property [7]. 

2. Framework 

Let the total system consist of a small system S and a large reservoir B. The 
Hamiltonian of the total system reads 

H = Ho + XHsB = Hs + Hb + XHsB, (4) 

where H^, Hb, and H^b are the Hamiltonians of system S, reservoir B, and the 
interaction between them, respectively, and A is the coupling constant. The cor- 
responding Liouvillians are defined as the commutators with the Hamiltonians, 

C = -i[H, ■]=Co + XCsB = Cs+J^B + A£sB. (5) 

Clearly, 

[Cs,Cb]=0. (6) 

We assume that the system Hamiltonian H^ admits a pure point spectrum, and 
the system Liouvillian £s is resolved in terms of its eigenprojections Qm, 

m m 

In this article, we discuss the derivation of master equations via Nakajima- 
Zwanzig's projection method, to implement the trace over the reservoir (^. We 
define two projection operators: V defined by Q and 

Q = l-V. (8) 

They are projection operators since they satisfy 

V + Q = l, V^=V, Q'^ = Q, VQ=QV = 0, (9) 

assuming the normalization of the reference state trB ^b = 1- One of the main 
issues of the present article is to discuss the choice of the reference state CIb ■ We 
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do not know at this moment which state should be taken as the reference state Ob- 
We only assume its stationarity here, jCb^b = 0, and will specify other features 
later. 

Notice that 

[P,/:s]=0, [Q,£s]=0, VCB = CBV = f), QCb = CbQ = Cb, (10) 

and that we can always accomplish 

VCsbV = 0, (11) 

by redefining the Hamiltonians H^ and H^b- 

3. Projection Method 

Let us now start the derivation of the master equation for psit) defined by (^. 
Consider the initial-value problem for the total system, 

j^p{t) = Cp{t), p{0)=po, (12) 

where the initial density operator po is not assumed to be factorized like Q- By 
projecting the Liouville equation (|12|) onto the two subspaces defined by V and Q, 
and noting the properties l\W\i and (fTT|) . one gets 

-Vp = CsTp + XVCsB Qp, (13a) 

at 

-Qp = C'oQp + \qCsbTp, (13b) 

respectively, where 

C'o = Co + XQCsbQ. (14) 

By formally integrating out the second equation and plugging the result into the 
first, one gets the following equation for the 7^-projected operator in the interaction 
picture [6], 



dt' 



-e-'^^^Vp{t) = X^ [ dt' e-^^^VCsBe^'o^^'^'^CsBVpit') + Xe-^^^VCsBe^'^^^Qpo- 
' Jo 

(15) 

This is an exact and non-Markovian master equation since the first term on the 
right-hand side contains a memory integral. The last term, on the other hand, 
represents the contribution arising from a possible initial correlation between sys- 
tem S and reservoir B. In fact, this term is absent when the initial state po is 
factorized like ^ and the reference state of the projector is accordingly chosen as 
^B = Pb • We are interested in the fate of this inhomogeneous term in the following 
arguments, when there is an initial correlation between system S and reservoir B. 
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4. Van Hove's Limit 

Let us make the Markovian approximation here, since the exact master equa- 
tion (|15j) is usually too complicated to handle for general systems and the Marko- 
vian approximation is valid for various practical situations. In this article, we 
apply van Hove's limit [8] for this purpose. That is, we take the weak-coupling 
limit A ^ and the long-time limit t — s- oo at the same time, keeping the scaled 
time T = X^t finite. This rescaling in time eliminates the memory, and the system 
exhibits a Markovian dynamics in the scaled time r. 

Now consider the density operator 

>p(r/A2), (16) 
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+ l.e~Csr/x^VCsBe'''o-/''Qpo, (17) 



with the initial condition 



p\->iO)=Vpo. (18) 

By integrating (fT7|) . one gets 

Jo Jo 

+ \r dr' e-^^-'/A V£sBe^o-'/A^ q^^, (19) 
"^ Jo 

which is arranged into the form [7] 

p!'^(r) =Vpo + Y. rdr'e^(— ")-'/A^;cW (^ _ r')p^^\r') +TW(r), (20) 

m,n ^ 

where 

K-^^ni-r) = VQn.CsBngHT)CsBQnV, (21a) 

TW(r) = X^VQmCsBni^\r)Qpo. (21b) 

m 

The important feature is that the kernel }Cmn{T) and the contribution of the initial 
correlation I^ > (r) are both governed by the operator 

f-r/X^ 

7^W(^)= / ^^Q^{C'o+iuJm)t_ (22) 



/o 
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Therefore, the van Hove hmit of this key operator TZm (t) completes the van Hove 
Umit of the master equation H2U|) . 

5. Main Theorem 

As we wiU see later, the van Hove limit of TZm (t) results in 

'J^i^Hr)^=^- .^ (r>0), (23) 

and it immediately leads us to the conclusion that, in van Hove's limit, the memory 
kernel lQnn{T) in (|21aj) is reduced to a Markovian generator 

^^^n{r) ^^ /CW = -VQ^C^B-^—^ ^CsbQuV, (24) 

while the correlation term X''^' (r) in ()21b|) disappears 

j(A)(^)^^0, (25) 

so that the reduced dynamics (^0]) becomes 

pW(^) ^^Ppo + V r dT'6^nlC^±pi{T')=Vp,+ r dr'Kpiir'). (26) 



The master equation in van Hove's limit reads therefore 



(27a) 
with the initial condition 

pM = Vpo = iiB{po](^^B- (27b) 

That is, even if the initial state po is not in a factorized form, but rather there 
is entanglement, or simply a classical correlation, between system S and reservoir 
B, all correlations disappear in van Hove's limit and system S behaves as if the 
total system started from the factorized initial state ()27b|) with a reservoir state 
Qb specified below. 

The key formula (|23j) . and therefore the master equation (|27j) in van Hove's 
limit, are proved under the following assumptions: 

(i) the initial (correlated) state of the total system, po, belongs to a single sector. 
That is, pq is given in the form 

Pq = Y,U{Is®VLb)lI (28) 
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where Is (8) ^b is the state which specifies the sector and Lj's are bounded 
operators; 

(ii) the state Ob is mixing, and therefore the reservoir Liouvilhan £b bears a 
simple eigenvalue 0. We assume in addition that the remaining part of the 
spectrum of Cb is absolutely continuous; 

(iii) the projection operator V adopted in (^5]) and (^Tj) must be defined in terms 
of r^B that is mixing and "contained" in the initial state po of the form (|2^ 



6. Mixing 

The state Q,b is called mixing with respect to the reservoir dynamics e^^* iff 

trB{X(i)yQB} = trB{Xe^Btyf^g} hl^ trBl^f^B} trBlrOB} (29) 

for any bounded superoperators X and Y of the reservoir B, where X{t) = 
e~^*Xe ^3* [9], One of the typical mixing states is the thermal equilibrium state 
of bosons at finite temperature Qb ~ e~^^^ [7], which often appears in the lit- 
erature. Among other interesting mixing states, there are nonequilibrium steady 
states, where system B consists of two reservoirs with different temperatures and 
steady current flows between them. Such a situation is also within our scope. 

Note here that a state pQ that refers [in the sense of (|28|) ] to two (or more) 
different mixing states of an infinite reservoir cannot be a physical state, since it 
is the superposition of states belonging to different inequivalent sectors. (Imagine, 
for example, the superposition of states with different temperatures.) Hypothesis 
(i) is therefore reasonable from this point of view. 

The important consequence of mixing is that it implies the existence of the 
point spectrum of Cb [9]. Let us consider the time evolution, driven by the 
reservoir free Liouvilhan Cb, of the expectation value of any operator D = ^ • Aj® 
Xj, starting from the initial state pQ of the kind ()28() . where Aj^s are operators of 
system S and ^j's bounded operators of reservoir B. Then, the mixing property 
(jk|) of Qb yields 

tr{Z?e^sVo} = J] J^trB[X,e^Bttrs{AjLi(ls f^B)^}] 

i 3 

^^^^ X]I^trB{Xjl]B}tr{AjL,(ls ® f)B)4} = tr[Z?(trBM ® Ob)], (30) 
« i 

and in this sense, we have 

e^^Vo ^=^ trBJpo} ® ^B (31) 

for any state po belonging to the sector of Is ® ^b- This shows that Cb bears one 
and the only simple eigenvalue (in the relevant sector), which is 0. It is essential 
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to note that, if Cb admits more simple eigenvalues, such as for an ergodic state 
r^B [9], it is not possible to derive a master equation as we will see later. Recall 
further that mixing does not exclude a singular continuous spectrum, as far as it 
is transient, i.e., Riemann-Lebesgue's lemma holds [9]. Hypothesis (ii) therefore 
assumes stronger conditions than mixing. 

Finally, Eq. (|cill) also shows that the eigenprojection Ho belonging to the simple 
eigenvalue 0, 

£Bno = no£B = 0, (32) 

acts on the state p in the relevant sector as 

Uop = tTB{p}®nB. (33) 

The projector V suggested in Hypothesis (iii) is nothing but this eigenprojection 
Hq: the proper choice of the projector is 

T' = Ho, (34) 

where the reference state Qb should be endowed with the mixing property and 
contained in the initial state as (|28)1. A wrong projector yields secular terms and 
prevents us from deriving a master equation, as we will see in the next section. 

7. The Key Formula 

Let us finally prove the key formula (jSHJ with the projector (|34|) . It requires a 
few steps, since the spectral property of Cq given in H14() is not known. We start 
by noting that Qe^o* = Qe^* Q and hence ||Qe^o*|| = ||Qe^*Q|| < ||Qp • ||e'''*|| = 1. 
One therefore obtains 

A"7^W(r) = A" / dt Qe(^o+*^™)« ^ o as A ^ (n > 2) (35) 

Jo 

in van Hove's limit, irrespective of the spectrum of Cq. 
Second, the convolution 

/ dt'e^°(*-*')Q£sBQe^o*' (36) 

Jo 

is bounded for any t, provided the projector is properly chosen as V = Ilo and the 
point spectrum of Cq is removed by the projection Q = 1 — Hq. See Ref. [7]. 
Now, by using 



e^o* = e^o* + A /* dt' e^°(*-*') Q£sB Qe^o*' 
Jo 



(37) 
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we expand the relevant quantity ((^ as 

j-t/X^ [-t/X^ ft 

7^W(r)=/ dt Qe(^«+^'^'")* + A / dt dt' e^^''+''^"^^^e-^°^' QCsBQe^'^'^' . 
Jo Jo Jo 

(38) 

Let us split the first term into two contributions, the discrete and the continuous 

spectra of Cq, by the projection operators XIq and 1 — Ho. For r > 0, one gets 

r/A2 

dt ge(^o+*'^m)i 



/A2 /.r/A2 



Jo Jo 



'0 ^0 



i(Wm-Uln)T / X^ _ -^ ^{Co+iuJm)T / X^ _ ^ 



;Q^0Qm + V QJloQn 7 ^ + Q(l - Ho)- 



(39) 

which shows that the only possible divergence of the relevant operator in ()22() 
in van Hove's limit A ^ comes from the point spectrum, of Co (i.e. the first 
term of the last expression). If V = Hq, however, the divergent term disappears 
due to QHq = 0: the contributions of the point spectrum are properly removed by 
Q = 1 — Ho. The correct projection operator 7^ = ITo allows us to arrange the 
second term in (|5H|l into 

rr/X'^ 

A 



r dt [ (it'e(^«+^'^™)*e-^«*'Q£sBQe^o* 
Jo Jo 



-^ lUJ^TlX^ I 





JLq + lUJm - 0+ Jq 



- ^ r , -^ ^^SB r dt 26(^0+^"-)*, (40) 

Co + lUJm - 0+ Jq 

and we arrive at the recurrence formula 



^{Co+iLiJm)r / X'-^ _ -|^ 



Co + ioJm - 0+ 

o rl^^ 

+ A- ^e*"™"/^' / dte^°("/^'-*)Q£sBQe^o* 

£o + 1''^m - 0+ Jo 

-^7^—7^ ^'CsB^L'^(r) (r>0). (41) 

The integral in the second term is the convolution H36() . which is bounded for 
r/A^ -^ oo since the point spectrum of Co is removed by the projection Q = 1 — Ho. 
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Therefore, by iterating the expansion (|41|) twice and by noting 

JCo+iuJm)t 

together with (|35jl . we derive the key formula (|23|) . 

In summary, the correct projector "P = Ho helps us to manage the subtle point 
spectrum of £b- With a wrong projector, we would fail to appropriately remove 
it, Eq. ()23() would not hold anymore, and TZm (t) would diverge in van Hove's 
limit. This is because the wrong operator projects the reservoir onto a wrong 
state, so that the system evolution is frustrated, and develops secular terms in 
r/A^ [like in Eq. P?])]. 

As stressed in Sec. 6, the reference state of the reservoir must be mixing in 
order that no discrete eigenvalue different from exists. Otherwise, the point 
spectrum (except 0) is out of control, since we do not know the eigenprojections 
other than Ho, and again the emergence of secular terms is inevitable. Moreover, 
it should be noticed that the mixing property of the reservoir is crucial even for a 
factorized initial state like ((2)). 

8. Concluding Remarks 

In this article, we described the derivation of master equations via Nakajima- 
Zwanzig's projection operator method in the presence of the initial correlations 
with a reservoir. In van Hove's limit, the initial correlation disappears and the 
total system behaves as if it started from a factorized initial condition. It was 
demonstrated that (even for a factorized initial condition) the choice of the pro- 
jection operator is crucial and a wrong projection provokes secular terms. The 
mixing property of the reservoir is essential to the derivation of master equations. 

The factorization of the initial state in van Hove's limit implies that the total 
system is factorized at any moment in the scaled time r, since the "initial time" 
i = has no particular status. It is actually possible to prove it, i.e. [7] 

lim Qpir/X^) = 0. (43) 

It is worth noting that the interaction between system S and reservoir B is not 
essential to the factorization: free evolution eliminates the correlation, and the 
reservoir relaxes into the mixing state ^b. Indeed, for any state po of the total 
system of the type H28() with a mixing state J7b, we have 

g£ot^p t^^ e^stnopo = e^s* trBJpo} ® ^B, (44) 

where the contribution of the continuous spectrum decays out due to the mixing 
property. [See Eq. (|3T|) .] In the rescaled time r, the factorization is very rapid 
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(order A^ in the scaled time r), and the total system looks factorized at any 
moment. The timescales of the factorization and the mixing are related not only 
to the spread of the initial correlation and the perturbation from the mixing state 
but also to the size of the relevant reservoir observables. These issues are discussed 
in detail in Ref. [7]. 
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